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Abstract. We study the elementary excitations of a transversely confined Bose-Einstein condensate in 
presence of a weak axial random potential. We determine the localization length (i) in the hydrodynamical 
low energy regime, for a domain of linear densities ranging from the Tonks-Girardeau to the transverse 
Thomas-Fermi regime, in the case of a white noise potential and (ii) for all the range of energies, in 
the "one- dimensional mean field regime", in the case where the randomness is induced by a series of 
randomly placed point-like impurities. We discuss our results in view of recent experiments in elongated 
BEC systems. 

PACS. 03.75.Kk Dynamic properties of condensates; collective and hydrodynamic excitations, superfiuid 
fiow - 05.60.Gg Quantum transport 



1 Introduction 

The rapid developments of coherent atom manipulation 
which has recently allowed to study atomic interferom- 
etry of Bose-Einstein condensate (BEC) on a chip |1I2| 
opens up the prospect of considering a whole set of new 
transport phenomena in BEC systems. This can be con- 
sidered as a new domain for studying the concepts issued 
form mesoscopic physics. As for the clean 2D electronic 
devices considered in this latter field, the BECs are gen- 
uinely phase coherent. Moreover, whereas interactions are 
difficult to model in mesoscopic physics, their effects in 
BEC systems are rather well understood and are expected 
to lead to a whole body of interesting phenomena: atom 
blockade [S] , perfect solitonic-like transmission over a bar- 
rier 131) non linear resonant transport breakdown and 
revival of Bloch oscillations [H], to mention just a few ex- 
amples. 

Coherent transport phenomena are of special interest 
in presence of disorder. Interference effects have then a 
prominent role, resulting, in the non interacting case, in 
weak or strong localization, as observed in many differ- 
ent fields (electronic or atomic physics, acoustics or elec- 
tromagnetism) . The influence of interaction on this phe- 
nomenon are of great interest (see, e.g., the review 7 ) and 
have recently been addressed in the case of repulsive two 
body effective interaction for BEC systems in Refs. |8l9j . 
In these latter two references, interaction effects have been 
shown to lead to genuinely non-linear phenomena that 
profoundly alter the usual picture of Anderson localiza- 
tion. 

In the present work, we also consider the influence of 
interaction on Anderson localization, but remaining at a 



linear level, by studying the propagation of elementary 
excitations in a disordered BEC system. These are small 
deformations of a static background and they can be -at 
leading order- described in a linear framework (neglect- 
ing phenomena such as Beliaev damping) . Interaction has 
nonetheless a prominent effect on the spectrum of elemen- 
tary excitations, which is phonon-like at small energy and 
becomes similar to the one of non interacting particles at 
high energy. The crossover between these two regimes oc- 
curs at an energy huj of order of the chemical potential fi 
of the system. 
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Figure 1: Lioc as a function of the energy huj of an ele- 
mentary excitation in logarithmic scale (fi is the chemical 
potential of the system, and is the value of L\oc when 
huj = fi). The curve has been drawn within the model 
used in Section^ employing formulas i)4(J)) and i)6'l|) . This 
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(the meaning of the parameters X and rtimp in this for- 
mula is explained in Sectional)- 
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Accordingly, the localization length Lioc of the ele- 
mentary excitations (i.e., the typical extend of a localized 
mode, see Section 2 below) is expected to be similar to 
the one of phonons at low energy {huj <C fJ.), and to the 
one of non interacting particles at high energy {fku » fi) 
|10| . The localization length of non interacting particles 
scales linearly with the energy (at high enough energy, 
see, e.g., whereas phonons in ID disordered system 

have a localization length which diverges as uj~'^ at small 
w, as typically observed in models of disordered harmonic 
chains |12j . in random layered media 13 , or in contin- 
uous models with random elastic properties ^1]. Hence, 
the localization length Lioc of the elementary excitations 
has the behavior illustrated in Figure 1, with a minimum 
at fiw ~ /i. The main purpose of the present work is to 
explicitly derive this type of behavior within several ap- 
proximation schemes and different models of disorder. 

The paper is organized as follows. In Section El we 
briefly present the model and the parameter range in which 
we are working, together with the Bogoliubov-de Gennes 
equations governing the dynamics of the elementary ex- 
citations. In Section Owe consider the large wave-length 
limit within an hydrodynamical approach. We consider 
a Gaussian white noise potential and show in particu- 
lar that in this domain, one obtains a u!~'^ behavior of 
Lioc. In Section 01 we consider an other type of disorder 
(randomly placed delta impurities) and work within the 
transfer matrix approach. In this regime we are able to 
work for all the range of energies and obtain an analytic 
expression for Lioc in the scarce impurities limit. This ex- 
pression matches at low energy the one obtained in Section 
within the hydrodynamical approach. Very interesting 
recent experiments have addressed the issue of transport 
in a disordered BEG jl5ll()ll7j and in Section Owe dis- 
cuss the relevance of our approach for analyzing some of 
the experimental results. Finally, some technical points 
are given in the Appendices. Appendix is devoted to 
the derivation of a formula allowing to determine the den- 
sity of state within the "phase formalism" employed in 
Section O In Appendix ^ we compute the transmission 
coefficient of an elementary excitation of energy fko over 
a single delta-like impurity. 

2 The model 

In this Section we present the basic equations describ- 
ing the elementary excitations of a one dimensional (ID) 
Bose-Einstein condensed gas in presence of disorder. The 
condensate is formed by atoms of mass m which interact 
via a two-body potential characterized by its 3D s-wave 
scattering length a > 0. The gas is confined to one di- 
mension by a transverse parabolic potential of frequency 
u!j_ and "oscillator length" aj_ = {h/muj^y/^ . There is 
no confinement in the axial (x) direction, but disorder is 
induced along the axis of the guide through a random po- 
tential U (x) whose properties will be specified in the next 
sections. 

In this Section (and also in Section 0} we restrict our- 
selves to the "ID mean field regime" [Ej corresponding 



to a density range such that 

(a/ax)^ < niDO < 1 , (1) 

where niD denotes a typical order of magnitude of the ID 
density n{x, t) of the system. The first of the inequalities 
Q ensures that the system does not get in the Tonks- 
Girardeau limit and the second that the transverse wave 
function is the ground state of the linear transverse Hamil- 
tonian, see, e.g., the discussion in Refs. |18I19| . We address 
the low density case (Tonks-Girardeau limit) and the high 
density case (transverse Thomas- Fermi) in Section 

In the ID mean field regime, the field operator is a 
function ^{x,t) which can be decomposed in the usual 
Bogoliubov way in c- number (the superfluid order param- 
eter) plus small terms describing the contribution of the 
elementary oscillations (see, e.g., Ref. [201: chap. 5). For a 
stationary condensate, the order parameter is of the form 
ipix) exp{— i/ii/fi} where ipi^) is real, and the Bogoliubov 
decomposition reads 

>F{x,t) = e-'f"'/^ ^^ipix) + 

5]K(x)6.e-'-'*+<(x)6te'-'i}, (2) 

where hy and 6^ are, respectively, the annihilation and 
creation operator of the v\h elementary excitation. In the 
following, we drop the subscript v for legibility. The order 
parameter verifies the Gross-Pitaevskii equation 

+ \u{x)+g^^i?{x)^^,{x) = yi^,{x) , (3) 

with (7iD = 1f\u}\_a [2112 2 23 . The functions u{x) and v{x) 
are solutions of the Bogoliubov-de Gennes equations (see, 
e.g., Ref. 1201, chap. 5) 

where 

ii = -^^^V{x) + 2g^y,^\x)-y.. (5) 
Irti dx^ 

In presence of a single elementary excitation of pulsation 
LO the density reads n{x,t) = \tp{x)\'^ + Sn{x,t) where the 
density oscillation is, at leading order: 

Sn{x, t) = il}{x) [u{x) + v{x)] e"'"* + c.c. , (6) 

where "c.c." stands for "complex conjugate" . In Section 
Owe use the notation 5n{x) for the quantity il)(x)[u{x) + 
v{x)\. 

In the absence of potential [/, the order parameter is a 

1/2 

constant %l^{x) = with fi — gyo no, the speed of sound 
in the system is cq = {ii/mY^^ and the healing length is 
C = fi/imco). 

Disorder is induced along the axis x of the guide through 
the random potential U{x). Denoting C/typ the typical 
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value of |C/(x)|, we work in the limit C/typ <C M- This regime 
is easily reached experimentally j 151 161 and is very rele- 
vant for our purpose because it corresponds to a range 
of parameters where Anderson localization is not blurred 
by effects connected to "fragmentation of the condensate" 

El- 

In a ID disordered system the excitations are expected 
to be localized around a point with an envelop decreasing 
exponentially with the distance to this point. This cor- 
responds to functions u, v and Sn behaving as exp{±7a::} 
when |a;| — > oo. 7 is a function of lo known as the Lyapunov 
exponent; it characterizes the localization properties of the 
system. Its inverse Lioc = is the localization length 
We determine the Lyapunov exponent of the system 
in Section 13 in the hydrodynamical regime huj <C ^. 

In Section 0] we approach the problem in a different 
-but equivalent- manner. The disordered potential is as- 
sumed to be non zero only in a finite region of space, 
between a; = and L. We consider an elementary exci- 
tation of pulsation lo incident on the random potential. 
The corresponding transmission coefhcient T through the 
disordered region is related to the Lyapunov exponent via 
7 — — i lim^^oo biT 11 . This is simply connected to 
the fact that the incident wave function decreases expo- 
nentially -at a rate 7- in the disordered region, and this 
corresponds finally to a transmission probability which is 
(within logarithmic accuracy) T ~ exp(— 27L). 

We note here important features of the localization 
properties of the elementary excitations. First, Eq. Q) ad- 
mits a zero energy solution for u{x) — —v*{x) = ip{x). 
Thus, whatever the disordered potential U{x), the excita- 
tion at w = is delocalized since tpix) extends to infin- 
ity. This implies that Lioc diverges as lo ~^ 0. Secondly, 
at w ^ 00 the high energy part of the spectrum is well 
described by a single particle description obtained by ne- 
glecting the coupling between the positive (u) and neg- 
ative (u) frequency components of the excitations (see, 
e.g., Ref. [20 , chap. 12). In this limit one can set u = in 
Eq. and the system is described by the Schrodinger- 
like Hamiltonian H ^ which localization length behaves 
as iioc c>c oj at high energy. Thus, as already anticipated 
in the introduction, we expect a behavior of Lioc similar 
to what has been drawn in Fig. 1. 



Thomas-Fermi approximation j25| leading to 



lfi-U{x) 



9lB 



(7) 



By reintroducing this ansatz in Eq. (PJ, one can easily 
show (provided C/typ is small compared to /i) that the 
Thomas- Fermi result Q is valid in the limit ^ <C 
where Tc fixes the length scale of typical variations of U 
(for instance this is the correlation length of the random 
potential). If besides, one considers the limit ^ <^ cq/lo, 
the density oscillations dn{x) obey the hydrodynamical 
equation j26l27| 



-uj'^Sn{x) 



dx 



c^(x)—Sn(x) 
dx 



(8) 



where c{x) = {[/i — U{x)]/m}^/^ is a local sound velocity. 

Disorder is induced along the axis of the guide through 
the random potential U (x) which is assumed to have zero 
mean. The case (U) 7^ can be treated with a trivial ex- 
tension of the present approach which is explained at the 
end of the Section. In the following of this Section, U will 
be approximated by a Gaussian white noise. The hypoth- 
esis of white noise is only valid if the wave length of the 
excitations is large compared to the correlation length 
of the true U (which is not a perfect white noise if we want 
the Thomas Fermi approximation to hold). Hence, in 
the present Section, we make the consistent hypothesis 
that 

2 TT Co 

C < < . (9) 

LO 

When the inequality @ is verified. Equations O and © 
are both valid and furthermore the approximation of the 
random potential by a white noise is sound. In the follow- 
ing we thus write 



{U{x)Um^[-] D5{x). 



(10) 



We now evaluate the localization length corresponding to 
Eq. ((HJ by means of the phase formalism (see Ref. jll|L 
We consider a real solution of (jSJ and define the functions 
a{x) and (3{x) by 



3 Hydrodynamical approach: hu <^ /j, 

The results obtained in this Section are derived within the 
ID mean field regime (0. As explained at the end of the 
Section, they can be easily generalized in the transverse 
Thomas-Fermi regime and even in the Tonks- Girardeau 
limit. 

In the present Section we only consider the low fre- 
quency excitations {fko <C /x). These involve large wave 
lengths (which are of order 27rco/w, when w — > 0) and 
accordingly, features at small length scale are not rele- 
vant in the potential seen by the excitations. In particular, 
the ground state order parameter can be evaluated in the 



a{x) 



Sn{x) 
6n* 



c^(x) da 



Co oj dx 



-jz- (11) 



In Hll|) the quantity 5n* is a typical value of Sn{x) which 
is introduced for dimensional purpose, but plays no role 
in the following [since Eq. is linear]. The functions a 
and P satisfy the following system of equations: 



da 
da; 



Co 



1 + ri{x) (3{x 



d/3 
da; 



Co 



.a(a;). (12) 



In the first of equations (|12|) . the term ri{x) is equal to 
U{x)/[iJL — U{x)]. In all the following we assume that C/typ 
is much smaller than /x, and we write ri(x) ~ U{x)/ fi |28j . 
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It is convenient to parametrize the functions a and f3 
in the form 

a{x) = r{x) cos6{x) , P{x) — r{x) s'm6{x) . (13) 

The functions 9{x) and r{x) describe respectively the phase 
and the envelope of the density oscillations Sn(x) [and ac- 
cordingly of u(x) and of In particular, the Lyapunov 
exponent is defined by 



Simple algebra allows to express the average {z)st = 
Im. ^st(-^) "^-^ under the following form: 



{z)st = Ju. 





12t2 + i6- 


/ exp 


6 ujS 



(22) 



7(w) 



lim 

X — ^oo 



(lnr(a;)) 



(14) 



It is convenient to introduce the quantity z — a/(3 be- 
cause, owing to the equality 

2lij 

In (x)^— z{x') dx' + In f3{0) - In sin^ 0{x) , (15) 
Co Jo 

and to the fact that the probability density of sin 9 (and 
thus also that of z = 0016*) becomes stationary (i.e., x 
independent) at large x |llj . one can write 



7 



— lim X 

Cq X — ^oo 



z{x')) dx' = - (2),t , (16) 

Co 



where (-z)st is the mean value of z in the stationary regime. 
This quantity is determined as follows. From (|12|l one sees 
that z verifies the following stochastic differential equation 



-^^ = l + z^ 
LU dx 



U{x) 



(17) 



Let P{z; x)dz be the probability that z{x) lies in the in- 
terval z, z + dz. From IjlTI) and (|10|l P verifies the Fokker- 
Planck equation (see, e.g., |Tll29 |) 



dP _ Lo d 
dx Co dz 



ujSdP 



(18) 



where 5 = D /co . The stationary regime corresponds to 
the case where dxP = 0. In this case, writing P = Pst{z), 
Eq. (dHJ yields 



{l + z'')P,t + 



LoSdP^t 

2 dz 



(19) 



where J^^ is an integration constant. The solution of H19|l 
is 



Pst{.z) = 



2J^ 

LU 5 



+ 00 

dt 



-{I + z^)t + zt^ 



(20) 



The value of Joj is fixed by the normalization of Pgt- One 
obtains 



J- 







/ exp 


6ujS 


m 



dt . 



(21) 



We are primarily interested in this Section in the small 
frequency evaluation of the Lyapunov exponent, because 
Eq. © is expected to describe the elementary excitations 
only in the domain hu/ <g; 1 . An expansion of the in- 
tegrals H21|l and H22|l in the limit ujS ^ yields, after 
reinserting in (|16|l : 



^ 15 /'ujS 
^ 16 \Y 



(23) 



Although the high frequency limit is not expected to be 
relevant in the hydrodynamical regime, we note for com- 
pleteness that when coS — > cxd one obtains 



Co 



1 - 



p 



V6 
uj5 



\ 2/3 



1/3 



2[r(j 



3/^(1) 



Ve 



4/3 



(24) 



The exact value of 7 -as determined numerically from 
Eqs. (|2U and is represented in Fig. 2 (solid 

line). 



0.5 - 




Figure 2: 7 as a function of lu in rescaled units. The 
solid line is the numerical evaluation of 7 using formulas 
nib]) and y22\) . The dashed lines are the small and large 
ujS approximations [Eqs. and y24\)]. 

The quantity Jui is also of interest for itself, because it 
gives informations on the density of states of the excita- 
tions. It is show in Appendix IXI that . if N{ijj) denotes the 
integrated density of state per unit length, one has 



Co 



(25) 
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From H21(l one gets the following expansions: 



N{uj) 



3a; r(|) 




when ivS 3> 1, and 



TTCo 



1 



5 

32 



lij5 

T 



(26) 



(27) 



when 0)5 <C 1. In the relevant regime of low excitation 
energies, the leading order in H27|) coincides the result in 
absence of disorder, where one has a linear dispersion re- 
lation ljj = colgl in the hydrodynamical regime. This con- 
firms what could have been already anticipated from the 
fact that 7 — > when w — > 0: the low lying excitations are 
poorly affected by the presence of disorder (the relevant 
small parameter being lij5). In particular, there is no trap- 
ping of the elementary excitations by the disorder and no 
Lifshitz tail in the density of state. This is linked to the 
fact that = constitutes what is called a "stable gen- 
uine boundary of the spectrum" in the book by Lifshits, 
Gredeskul and Pastur (see Ref. ^J, section 7.3). 

The results presented in this Section have been ob- 
tained for a random potential with zero mean. They are 
very easily adapted to the case {U) ^ 0: it suffices to 
write U{x) = ([/) + Ui{x), and to define = fi — (U), 
ci = (/ii/m)^/^, ^1 = h/mci, 6i — ^fD/ci. Then, all the 
results presented from Eq. Hll() to Eq. H27(l remain valid 
provided U{x), ^, cq, ^ and 5 are replaced by the similar 
quantities with subscript "1" , with the coefhcient D being 
now defined by {Ux{x)Ui{Q)) = {h^/mf D5{x) [instead of 

The present hydrodynamical approach is very interest- 
ing because it has natural extensions out of the ID mean 
field regime defined by Eq. (Q. For high linear densities, 
when niDfl 3> 1, one reaches the "transverse Thomas- 
Fermi regime" also named "3D cigar" in Ref. JHI- In 
this regime the system cannot be considered as truly uni- 
dimensional. However, the lowest branch of the spectrum 
corresponds to excitations that are isotropic in the trans- 
verse direction, and, as shown by Stringari in Ref. 30 , 
they can still be described within the hydrodynamical ap- 
proach. In this case, averaging the 3D hydrodynamical 
equation over the transverse direction, one gets a ID equa- 
tion of the form where the local sound velocity c{x) 
is now taken to be c{x) = — U{x)]/m\^/'^ . So, all 

the results presented from Eq. Hll|l to Eq. H27|l remain 
valid provided ^, cq and ^ and are replaced by /x' = /i/2, 
c[) = (/i7m)i/2 and = h/{mc'o). 

The low density regime mua <C {a/a±)'^ (Tonks-Gi- 
rardeau) can also be studied within the hydrodynamical 
framework (see for instance Ref. ^H])- In this case one has 
/i = {tt Huq)^ /2m, Co = (2/i/m)-'^/^ and Eq. (jSJ is replaced 

by 

—uj'^Sn{x) — — < c{x)— c{x) Sn{x) 



dx 



dx 



(28) 



with the local sound velocity being defined by c{x) — 
{[— [/i — U{x)]}^^^. In the present case we define [instead 
of'tni)] 



a{x) 



2{x) Sn{x) 



Co Sn* 



c{x) da 
dx 



Writing a = r smO and (3 = r cos 9, one obtains: 

dr ^ d^ uj 

dx ' dx c{x) 



(29) 



(30) 



From the second of these equations, in the limit where 
l/c(a;) = Cq ^[1 + ^U{x)/n], one can show that the phase 
0{x) has a Gaussian distribution of the form 



Q{0;x) 



1 



y^27rLL>^xS/co 



exp ■ 



CO ' 



2uP-xblcQ 



} , (31) 



with (5 = [W- l2m,[ifDlcQ. 

The first of Eqs. H30|l is more interesting. It shows that 
the envelope of function ol(x) remains exactly constant. 
Assuming that the localization properties of a[x) are the 
same than those of bn(x) .3I4, this equation points to the 
absence of exponential localization in the hydrodynamical 
limit of the Tonks-Girardeau regime. 



4 Transfer Matrix approach 

In this Section, we study Anderson localization of the ele- 
mentary excitations of a Bose-Einstein condensate with an 
other type of disordered potential and in a framework dif- 
ferent from the one used in the previous Section. Namely, 
we study the transmission through a disordered region of 
extend L, in the ID mean field regime (Q), by means of a 
transfer matrix approach for a disordered potential: 

U (x) = gimp ^ (5(a; - x„) , where 5imp = A /i <^ . (32) 



U(x) describes a series of static impurities with equal 
intensity and random positions Xn- The peak intensity is 
measured by the dimensionless parameter A. We consider 
here the repulsive case A > 0. The x^'s are uncorrelated 
and uniformly distributed with mean density riimp. In this 
case <U[x)>— gi„,pnimp and <f/(xi)C/(x2)> — <C/(xi)> 



X <\J{xi)>= {h^ lmYDb(xx~x-2),^\t\\D 



From what is known in the case of Schrodinger equation, 
this type of potential is typical insofar as localization prop- 
erties are concerned Besides, it has recently been pro- 
posed to implement a very similar type of random poten- 
tial by using two different atomic species in an optical 
lattice 122. 

The static background is deformed around each im- 
purity over a distance which is at most of order ^. We 
consider the regime where this deformation does not ex- 
tend to the nearest impurity (ni^pC 1 ESI)- In this 
case, the propagation of an elementary excitation in pres- 
ence of the disordered potential U(x) can be treated as 
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a sequence of scatterings over isolated perturbations. Be- 
sides -as shown in Appendix El" both the scattering of an 
elementary excitation over such a perturbation, and its 
propagation between two successive impurities (separated 
by a distance £) are, in this regime, described by a 2 x 2 
transfer matrix, denoted respectively Tx and To{i) with 
(see, e.g., [SI) 



described by the transfer matrix TJv which is the product 



' " ' -rx/tx 1/tx 



Me) 



_(ilti 



o%A„) • (33) 



rx and tx in Eq. (|33|l are the transmission and reflex- 
ion amplitudes of an elementary excitation with energy 
huj across the background deformation induced by a sin- 
gle delta-like impurity. Their dependence on A and uj is 
determined in Appendix IbI [Eos. ifTTjl and (|SHJ]- The scat- 
tering states we choose for writing the matrices Tx and 
7o(^) are the one introduced in this Appendix. They are 
pictured in Fig. 3. 




Figure 3: Scattering channels used for writing the trans- 
fer matrices Tx and Tq{£) of i)^^) . In the case of unper- 
turbed motion over a length £ one has Xout — a;in = ^- In the 
case of scattering by a delta peak located at (a;in -l-a;out)/2, 
one should take f <C a;out — xin ^ ^Tmp- 

The coefficients u^j and Vuj in Figure 3 are chosen in 
order to make the incoming and outgoing channels iden- 
tical to these appearing naturally in Appendix ^ when 
considering the scattering of an elementary excitation by 
a single impurity. One thus takes 



9? 
2 

2 



Co? 



(34) 



where q is defined in Eq. H5t)|) . In the case of scattering 
by an impurity, this corresponds indeed to the scattering 
channels defined by Eqs. (jHSJ, and l|K5|l . In the case 
of free motion over a length £, it is easy to see that these 
scattering channels correspond to a matrix Tq{£) such as 
defined in Eq. H33|l with 



to(^,w) = e*(«^-2") , where e" 



(35) 



Then, the scattering by a series of N delta peaks sepa- 
rated by distances £i = X2 — xi, £n-i = xn — xm-i, is 



Tn=Txx To{£n-i) X Tx... x To(^i) x Tx 
Tjv defined in Eq. H3t)|l is of the general form 



Tn = 



1/% ^n/^*n 

-TN/tN l/tN 



(36) 



(37) 



Eq. H37I) is used for computing the reflexion and transmis- 
sion amplitudes (rjv and tj^) of the elementary excitation 
over the potential H32|) . The transmission probability over 
this potential is Tjy — I^Arp. 

As discussed at the end of Section 2, the analogous 
of the Lyapunov exponent already computed in Section |51 
[Eq. (fTHl ] is here defined as 



7 = — lim 



''imp 



(ln|ijv|)-- lim 



' ''imp 

2N 



(38) 

We calculated 7 numerically, by a Monte Carlo averaging 
over 50 realizations of the disorder, taking N — 2000 2^. 
The result is shown in Figure 4 for A = 1 and rtimp S, = 
0.02. In the present model the lengths £i — Xi+i — Xi are 
independent, Poisson distributed, random variables with 
P{£) = nimp exp{— ^ Tiimp}- Thus, for a fraction of lengths 
equal to nimp £, the transfer matrix approach fails because 
the distance between two successive impurities is smaller 
than ^ |3f)j . This is the reason why we consider a rather 
small value of density of impurities: for the chosen value 
"impC = 0.02, only 2 % of the distances violate the crite- 
rion of applicability of the transfer matrix approach. 




^ 03/ 10. 

Figure 4: 7 as a function of lo in rescaled units. The 
plot is drawn for A = 1 and ^nimp = 0.02. The dots are 
the results of the numerical simulation and the solid line 
is the analytical result from Eq. il4(J)) . The inset displays 
a blowup of the figure at low energy. 



As shown in Ref. |38| . in the limit nimp <SC q, one 
can obtain an analytical estimate of 7. From the relation 

Tn+1 = Ta X To{£n) X Tn one gets 



(ln|tjv+i|) =ln|fA| + (ln|<Ar|) 

\n\l + rxr*N^-^tliiN)\ 



(39) 
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In is typically of order n-^^p, and in the limit riimp ^ g, 
one may assume that the phase of to(^Af,^) given in H35|l 
is uniformly distributed in [0, 27r]. Then, the last term of 
the r.h.s. averages out to zero |37I38| . This yields 

7 = -ni„pln|tA|=-^lnTA, (40) 

where we recall that the explicit expression of t\ is given 
in Eq. (|^ and T\ = I^aP- Formula iPHjl corresponds to 
the solid line in Fig. 4. The agreement with the result of 
the numerical simulation is very good, even at low en- 
ergy, as shown in the inset of the figure. This is not a 
surprise because the breakdown of H4U|I is expected only 
at extremely low energies for the present value njmp : when 
(7 <,ni,„p, i.e., hio/^ ~ ^ q <0.02. For larger values of nimp, 
the good agreement of Eq. H4U|) with the numerical data 
is limited to a smaller range of energies, mainly because 
the transfer matrix approach fails. 

From Eq. H60|) . in the limit of small oj and A <C 1 
formula (|^ yields 

(^-j r^i^p. (41) 

The precise range of validity of formula H41I) in the energy 
domain is expected to be ^Jiimp <C hto/fj, <^ 1; the first 
inequality ensures that (|40|l is valid and the second that 
H60(l is applicable. The accuracy of formula H41|l is tested 
in Fig. 5 in the case A = 0.2 and ^ ni^p = 0.03. As al- 
ready seen on Fig. 4, one notices on this Figure that the 
restriction ^ nimp huj/ ^ turns out to be of no practical 
importance. 



0.06 




Figure 5: 7 as a function of lo in rescaled units. The 
plot is drawn for A = 0.2 and nimp = 0.03. The dots are 
the results of the numerical simulation and the solid line 
is the analytical result from Eq. I^4(]\) . The dashed line is 
the approximate result H41}) . 

Formula (|41|l is interesting because it is identical to 
the first term of expansion (|23|l which has been obtained 
in Section|31in a completely different framework, and this 
permits to bridge the gap between the hydrodynamical 
approach and the present transfer matrix method. As just 
mentioned, formula H41I) is restricted to small values of A, 



but the approach of Section 13 is similarly limited to the 
domain [/typ ^ M- Also, Eq. H23|) is restricted to small val- 
ues of ujS. But in the present case ujS — ^nimp A^(?kj//i) is 
very small, even if /kj ~ /x, so the restriction ujS <^ 1 turns 
out to be of no practical importance here. Also, the results 
obtained in the present Section correspond to a potential 
with ([/) — A /i^ nimp 7^ 0. However, the comparison with 
the results of Section |31 is possible with the rule given at 
the end of this Section for treating the case of a potential 
with non zero mean. In this case the first term of expan- 
sion modifies to 7 = |A^nimp(fiw/^)^(l — A^nimp)"^- 
The correcting term (1 — A^nimp)"^, due to the non zero 
average of the potential, gives an undetectable modifica- 
tion of the result (the relative difference with (|41|l is of 
order 0.18 % in the case of Figure 5). 



5 Discussion and conclusion 

In this paper we have studied Anderson localization of 
elementary excitations in a ID BEC system. Emphasis 
has been put on the determination of the localization 
length which has been determined in Section |31 using the 
"phase formalism" in the hydrodynamical approach (valid 
for huj <C /i) and in Sectional using a transfer matrix ap- 
proach valid in the whole energy domain in the ID mean 
field regime (provided nimpf <SC 1). Results from the two 
approaches match within the appropriate limit. The hy- 
drodynamic approach has the advantage of being able to 
deal with a large range of linear densities, ranging from the 
low density Tonks-Girardeau regime to the high density 
transverse Thomas-Fermi regime. In particular the puz- 
zling absence of localization at low energy in the Tonks- 
Girardeau limit deserves further studies. 

Our findings can be tested in realistic experimental se- 
tups. Up to now, 3 experiments, lead at Firenze, Orsay 
and Hannover, have been done which all use similar con- 
figurations |15I16I17] . Each of these experiments involves 
an elongated cigar shaped condensate in a magnetic trap 
with an optical speckle pattern creating the disordered 
potential !39J . The experimental random potential has a 
non zero mean value, and the experiments are done in 
the transverse Thomas-Fermi regime. We can thus study 
localization in this configuration using (for excitations of 
energy small compared to the chemical potential /i) the 
above hydrodynamical approach of Section |31 adapted as 
explained at the end of this Section (replacing in all the 
formulas ^ by ^'^ = ^ - (U), c[ = {fj,[/m)^/^, etc.). 
One writes U{x) — {U) + Ui{x). The auto-correlation 
{Ui{x)Ui{0)) has a typical range Tc which is in all the 
cases much larger than the healing length ^: Tc = 20 fim 
and ^ — 0.35 fj,m in the Firenze experiment; Tc — 5.2 /im 
and ^ = 0.16 /xm in the Orsay experiment |40| : Tc ~ 7 /im 
and ^ = 0.3 /im for iV = 8 x lO'' atoms at Hannover. 
The condition @ is fulfilled provided the pulsation lj of 
the excitations is much lower than 2TTCi/rc (which, for in- 
stance is equal to 2tt x 340 Hz for {U)/fi — 0.2 in the 
Orsay experiment). In this regime, the potential can be 
approximated by a white noise with a coefficient D and a 
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correlation radius Tc defined by 

O D = J^{Ui{x)Ui{0)) dx ^ {Uf . (42) 

Following the procedure explained in Section 01 this leads 
to 

The experimental configuration which is closer to the one 
considered in the present paper is the one of the Firenze 
group jl5j which has studied elementary excitations of an 
elongated condensate in presence of a speckle pattern. The 
discrete excitation modes in elongated systems are sim- 
ilar to the continuous ones of infinite systems we have 
described in the present article only in the case of high 
quantum numbers (see, e.g., Ref. j^, chap. 12). Unfor- 
tunately, only the low lying dipole and quadrupole modes 
have been studied in Ref. ^H] • We nonetheless discuss this 
experiment using our results, keeping in mind that we can 
only provide rough orders of magnitude. 

The data of the Firenze group are presented in a way 
more easily analyzed within the model of random delta 
peaks of Section^ However, in the regime where Eq. H43|l 
is valid, all models are expected to yield the same re- 
sult, as argued in Section |3| and verified in Section ^ 
The only relevant parameter being the parameter D, or 
equivalently Tc [which is related to D by (|42|l ] . Within 
the model of random (5-peaks one has (U) = A/z^nimp 
and D = nimp(A/C)^ yielding = n^^p = 20 i^m 
The chemical potential in the Firenze experiment is /x = 1 
kHz and the excitations considered are the dipole (j^i — 
8.74 Hz) and quadrupole {1/2 = 13.8 Hz). For a disorder 
such that (C/)//i = 0.1 (which is typical in this experi- 
ment) the dipole excitation corresponds to a localization 
length Ljqp = 4.1 mm, whereas for the quadrupole one 
gets Ll^ = 1.6 mm {L^^ = ^L^^ since 1^2/1^1 = VW^)- 
We also note that higher excited modes having frequency 
i/n = ^\/n{n + 3) 41 30 have lower localization lengths: 
^\oc — n(n+z) ^loc- ^\oc bccomcs Comparable with the the 
typical axial size of the condensate (110 /im) for 71 ~ 10 

m 

A precise plot of the oscillations of a dipole mode is 
presented in Ref. in the case {U) j [i = 0.06 which cor- 
responds to a limit we can address using Eq. 
An experimental estimate of the value of the localization 
length can be obtained by fitting the experimental data 
with a sinusoidal oscillation at frequency v\ with a dam- 
ping exp{-2X(t)/L^^P}, where X{t) = ^Av^t is the dis- 
tance traveled by the dipole mode for an oscillation of 
maximal amplitude A. From the data presented in Ref. 
[TK] we obtain l!^^^ ~ 1.7 mm. This does not agree with 
the value L\^^ ~ 15 mm obtained from Eq. H43|l in the 
case (U) I [i = 0.06, but we recall that we do not expect 
the dipole mode to be equivalent to an excitation of an in- 
finite system. Thus, the damping observed in the Firenze 
experiment '15' cannot be accounted for by a model of 



infinitely long condensate with no axial trapping. Quan- 
titative theoretical description of this experiment should 
take the axial trapping fully into account. We nonetheless 
hope that the experimental study of higher excited modes 
could directly confirm the result (|43|l . 

It is also interesting to discuss the expected localiza- 
tion length in the Orsay experiment 16 , where the prop- 
erties of the random potential are well characterized. In 
this experiment, the potential is Poisson distributed with 
a mean value {U) which is a fraction of the chemical poten- 
tial (/i = 4.47 kHz). Taking {U) / ji — 0.2, and for instance 
u; = ojj, = 27r X 6.7 Hz (corresponding to the dipole excita- 
tion) one obtains L\oc — 11-8 mm. Besides, if one in able 
to generate excitations with w ~ 6 x w^, one still remains 
in the hydrodynamical regime and the above value of Lioc 
is decreased by a factor 36, becoming of the order of the 
axial size of the condensate (300 /im in the Orsay experi- 
ment TBI). We recall that the present approach does not 
strictly apply for low lying excitations of a trapped con- 
densate, but it is nevertheless interesting to get an esti- 
mation of the typical length scale for observing Anderson 
localization experimentally. 

We note that the Firenze 01] , Orsay TF and Han- 
nover (17) groups observed a saturation of the expansion 
of a condensate in a disordered potential. In the 3 ex- 
periments this phenomenon has been interpreted (see also 
45 46 ) as being due to the trapping of the wings of the 
condensate by the large peaks of the speckle potential, 
with no relation to Anderson localization. We hope that 
in the near future, new experiments will be able to di- 
rectly address Anderson localization of elementary excita- 
tions in transversely confined Bose-Einstein condensates, 
in configurations corresponding to the scenario analyzed 
in the present work. In this case, our study indicates that 
localization is more easily achieved for excitations of en- 
ergy of order ^ (see Fig. 1) created for instance through 
Bragg spectroscopy ^]. This range of energy is out of the 
hydrodynamical regime presented in Section|31 but the ap- 
proach of Section 0] allows to get a quantitative estimate 
of iioc in this case (for the ID mean field regime). 
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Paris XI et du CNRS, UMR 8626. 



A Appendix: Density of state within the 
phase formalism 

In this Appendix we briefly demonstrate Eq. (|25|) following 
a similar demonstration in Ref. TT . We first demonstrate 
that the phase 9 defined in (|13(l is a monotonic function of 
Lo. This can be shown by introducing the auxiliary variable 
y — ~cqz/uj. Expressing (|17|l in terms of the variable y, 
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differentiating with respect to uj and then integrating the 
resulting equation one gets 



dy 2(jj 
duj c; 



^ /dx'y2(:E')exp{— / y(x")Ax"} > 



.(44) 

Thus z = —Luy/co is a decreasing function of u, and 6 is an 
increasing function of lo (since is a continuous function 
and dz/d9 = —1 — z'^). verifies the equation 



.1 + sin^.^ 



(45) 



The fact that is an increasing function of u immedi- 
ately implies that the number of eigenmodes [solutions 
of © or equivalently of H12() ] with pulsation between 
and to, verifying the boundary condition cot(?(0) — Oq and 
cot 9{L) — 9l coincides with the number of pulsations 
w' G [OjCj] for which the accumulated phase 9{uj',L) as 
determined by H45(l with the initial condition 0{lu', 0) = 9o 
verifies 9{uj' , L) = 9^ + mn. This number is equal to 



E 



9{uj,L) - 61(0, L) 



E 



Oioj, L) - 9o 



(46) 



where E(a;) denotes the integer part of x. Passing to the 
limit L — > oo and allowing for the fact the number of 
states is a self averaging quantity one obtains 



N{uj) 



lim 



ttL 



(47) 



But the number (|46|l is also seen to be the number of 
times where the variable 9 equals zero modulo tt in the in- 
terval [0, L]. This stems from the fact that 9{x) can change 
interval [nn, (n + only toward a higher interval and 
cannot go backward to a lower interval, because, as seen 



from 103), d9/da 



uj/cq > 0. One may thus write 



N{uj) ^ lim — - 

L-+CXD Co 1^ 



dx Q"'''{0;x) , 



where 



^rcd 



(61; x) = (^S{9{x)-nn 

neZ 



(48) 



(49) 



is the probability density for the reduced phase. Owing to 
the fact that Q'^°'^ reaches a stationary (i.e., x indepen- 
dent) distribution Qlf{9), ^ yields 

N{cu) = - Q-'i(O) = - lim (1 + z^) P,,{z) , (50) 



Co 



Co 



where the last equality follows from the relation z — cot 9. 
The explicit expression (|20|l of Pgt evaluated at large z 
then yields the desired result (jSlJ. 



B Appendix: Transmission through a single 
delta peak 

In this Appendix we determine the transmission and re- 
flexion amplitude of an elementary excitation of energy 



h Lo incident from the left on a delta-like impurity located 
at x = 0. These coefficients have already been obtained 
by Kagan et al. in the case of a barrier of finite width 
|48j . In the present case the impurity interacts with the 
atoms forming the condensate via a potential A/x^(5(x) 
with A > 0. The condensate is deformed near the impu- 
rity and the order parameter reads 



ip{x) — tanh(|x/i^| + a) , with 



a — — sinh 
2 



(51) 



This form of ^{x) corresponds to two portions of black 
solitons matched together at a: = in order to satisfy 
the condition ^ [ip'{0~^) - V''(0~)] = 2 A V^(0). Far from the 
impurity (at x — > ±oo), the background is not perturbed 
and an elementary excitation of energy huj has a wave 
vector q such that uj = co(7(l-l-(7^^^/4)^/^, and is described 
by {u{x),v{x)) = exp{iqx){ui^,Vuj) where ~ by Eq. Q - 
the constants u^^ and are related by 



huj 
1 



= 



(52) 



The background is deformed near the impurity [as de- 
scribed by (|51|l ] , and in this region the form of the wave 
function of the elementary excitation is affected in a non 
trivial manner. However, one still has an analytical de- 
scription of the excitations around the stationary profile 
(I51|l because the expression of the excitation around a soli- 
ton is known (it is given by the squared Jost functions of 
the inverse problem 0^1) see also Appendix A of Ref. ISUjl. 
Thus one can write the appropriate incoming, transmitted 
and reflected modes of the problem. It is important how- 
ever to realize that the system has also evanescent modes 
localized around the impurity More specifically, the 
scattering process of an excitation of energy fko incident 
from — oo is described by 

S^-lix) = Ainc S;(-a;) + Arei ^,(-x) + ^^p{-x) , 

(53) 

when a; < 0, and 



(54) 



when a; > 0. The indexes "inc" , "ref" , "tra" and "eva" 
correspond respectively to incident, reflected, transmitted 
and evanescent channels. The expression of Sk{x) {k — q 
or ip) in and lfM|) is 



^fe(x) 



Akx 



M + _^+,tanh(| 



cok 



+ itanh(| -|- a) 



(55) 



and the quantities q and p are wave vectors related to uj 

by 



q^^V2[^/ihIJM' 

p^ = V2[^{f^Ju:r^ 



1-1 



1 + 1 



1/2 



1/2 



(56) 
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The wave functions defined in Eqs. H53|l and H54|l are 
the most general solutions of corresponding to an el- 
ementary excitation of energy fko incoming from the left 
and scattering on a potential U{x) = ^S{x). In partic- 
ular, the incident, transmitted and reflected components 
of (|53I54II all verify 1)5211 far from the impurity. The as- 
sumption nimp^ ^ 1 made in Section 01 ensures that the 
evanescent mode Sip does not reach the nearest impurity 
51 . This is the reason why the scattering on potential 
1)32(1 can be described via a transfer matrix approach us- 
ing only 2x2 matrices. 

The matching at a; = corresponds to S'^^^(O) = 
H(+)(0) and dS(+Vda;|p - dS'(-)/dx|„ = 2X^-^S{0). 
This yields a system of 4 linear equations determining the 

coefficients A^ci, Ara, ^cva and Actl in terms of ^inc- A 
tedious but straightforward computation yields 



tx = 



and 



where 



^ic 



2 + iq^ tanh(2a) A* 
-2 + ig$tanh(2a) ^ A 



2 -f ig^ tanh(2a) A* 
-2 + iq^ ta.nh{2a) A 



(57) 



(58) 



A 




(59) 



2?(p + iQ) tanha 



The transmission probability Tx = I^aP has the asymp- 
totic form Ta ~ 1 — \^ii/{2fvjj) when uj oo, and in the 
opposite small energy limit {huj <^ fi) one has 



A^O 



2jl 

Xhuj 

2n 



tanh a 



tanh(2a) 



(60) 



A typical behavior of Tx as a function of a; is plotted in Fig. 
6. The transmission probability is 1 at small frequency. 
This anomalous behavior of the transmission at small en- 
ergy has already been noticed in Ref. 0H] in the case of 
a barrier of finite extend. It is also in agreement with the 
findings of Ref. (50) where a dark soliton with velocity 
Vsoi Co (and thus reaching the limit where it becomes a 
mere density perturbation, i.e., a phonon, which is an ele- 
mentary excitation with q ^ 0) was shown to pass over an 
obstacle without radiating energy, i.e., without refiection. 



0.8 - 



0.6 - 



0.4 



0.2 



X = l 



X = 2 



1 2 . , 3 4 5 

Figure 6: Transmission probability Tx across a poten- 
tial U{x) — Xfi^S{x) as a function ofuj (in rescaled units). 
The solid line is the exact result )57)) in the cases A = 1 
and X = 2. The dashed line are the corresponding small X 
approximations i)6'l)) . 

The exact formula for Tx [from 1)57)) ] is compared on 
Figure 6 with an approximation valid for all lu when A ^ 
1: 



1 



1 



(61) 



It is seen on the Figure that this approximation is rea- 
sonably accurate already when A = 1. More precisely, for 
A = 1, the relative error due to the use of Eq. 1)61)) is lower 
than 3 % (the error is maximum around hio ~ 0-7 fi); for 
A = 0.5, this errors is lower than 0.5 %. 
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